Non-linear viscous saturation of r-modes 
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Abstract. Pulsar spin frequencies and their time evolution are an important source of information on compact stars and their 
internal composition. Oscillations of the star can reduce the rotational energy via the emission of gravitational waves. In 
particular unstable oscillation modes, like r-modes, are relevant since their amplitude becomes large and can lead to a fast 
spin-down of young stars if they are saturated by a non-linear saturation mechanism. We present a novel mechanism based on 
the pronounced large-amplitude enhancement of the bulk viscosity of dense matter. We show that the enhanced damping due 
to non-linear bulk viscosity can saturate r-modes of neutron stars at amplitudes appropriate for an efficient spin-down. 
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Compact stars are so dense that they could contain de- 
confined quark matter. To reveal their decomposition re- 
quires to connect microscopic properties to macroscopic 
observables. R-modes are oscillations modes of rotat- 
ing stars that are unstable under the emission of gravi- 
tational waves [1, 2, 3, 4]. In case of neutron stars they 
could spin-down young stars within a short time inter- 
val or limit the rotational rate of old stars that are spun 
up by accretion. However, since viscous damping can- 
not stop low amplitude r-modes in certain instability re- 
gions at high frequency they have to be saturated by some 
non-linear damping mechanism. We propose the strong 
suprathermal enhancement of the bulk viscosity at large 
amplitudes [5, 6] as a viable saturation mechanism and 
in contrast to previous studies [7] find that it can stop the 
growth of r-modes of neutron stars at amplitudes that are 
large enough for an efficient spin-down but small enough 
that viscous damping could dominate competing satura- 
tion mechanisms [8, 9]. 

The bulk viscosity of dense matter provides a mea- 
sure for the energy dissipation in a compression and rar- 
efaction cycle. It is maximal when the external oscil- 
lation frequency matches the time scale of the micro- 
scopic interactions that cause the dissipation, whereby 
the dominant interactions in the case of star oscillations 
are slow weak processes. The external density fluctuation 
An induces a corresponding oscillation of the difference 
of chemical potentials jUa which would vanish in weak 
equilibrium. Similarly, the rate of weak interactions van- 
ishes generally for fully degenerate matter in equilibrium 
and becomes finite either due to finite temperature effects 
or due to such deviations of the chemical potentials from 
their equilibrium value. The net rate for the considered 
weak process r<°' takes the general form 
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where the coefficients %j of the terms that are non-linear 
in the oscillating chemical potential difference and 
their number N are determined by the respective weak 
process. In this work we study neutron matter with an 
APR equation of state [10] at densities that are not high 
enough to allow direct Urea reactions, so that only the 
modified Urea process involving a bystander nucleon 
n + n^n + p + e~+v e , •••is kinematically allowed. For 
this process N = 3, corresponding to a strong non-linear 
dependence on T and fi& = jj.„— ji p — fi e . 

The bulk viscosity features three distinct characteris- 
tic regions [5, 6]. In the subthermal regime at low ampli- 
tudes, where /Ia <C T, jUa is linear in An and the viscos- 
ity is independent of the amplitude. The suprathermal 
regime ju^ > T is divided into a part where jUa is still 
linear in An but the viscosity strongly rises and another 
where the rise of jj,^ becomes weaker due to non-linear 
saturation effects and the viscosity decreases again. Pre- 
viously mainly the subthermal limit jUa <C T has been 
studied, but since the r-mode rises exponentially it even- 
tually reaches the suprathermal regime that will be stud- 
ied in the following. 

At sufficiently low temperature and large frequency 
the bulk viscosity allows a general analytic approxima- 
tion that covers both the subthermal and the part of the 
suprathermal regime where jUa is linear in An/h [5, 6, 7] 
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where C is a susceptibility describing the strongly inter- 
acting state of matter which in the case of neutron star 
matter takes the form 

= 4(l-2x)(4^-f) (3) 
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in terms of the symmetry energy S, the baryon density n 
and the proton fraction x. 




FIGURE 1. Left panel: The amplitude of the chemical potential oscillation H&/T as a function of the driving density oscillation 
amplitude An/n for different forms of dense matter. Right panel: Comparison of the bulk viscosity of the different forms of matter 
as a function of the density oscillation amplitude. In both plots, the frequency is a = 8.4 kHz, corresponding to an r-mode in a 
millisecond pulsar, and characteristic values T — 10 y K of the temperature and it = 2/jq for the density are chosen. The solid curves 
show the result for interacting matter, the dashed curved for non-interacting matter, whereas the dotted lines denote the analytic 
approximation eq. (2); for more details see [6]. 



The validity of the approximation eq. (2) can be seen 
from fig. 1, where on the left panel the amplitude of 
the chemical potential oscillation /J.&/T is plotted against 
the corresponding density amplitude An/n. For the dif- 
ferent forms of dense matter shown there the suprather- 
mal regime is reached for amplitudes An/n = (9(0.01). In 
particular, in the case of hadronic matter with modified 
Urea interactions the linear regime, where ji^ ~ An/n, 
extends at large frequencies basically over the entire 
physical range of density amplitudes An/n < 1. Corre- 
spondingly, the approximation eq. (1) shown by the dot- 
ted curve on the right panel of fig. 1 reproduces the full 
result shown by the solid curve favorably in this case. 

We consider in the following an exemplary 1.4 M Q 
neutron star model obtained as a solution of the TOV 
equations [11] using the APR equation of state [10]. The 
analysis of the bulk viscosity damping time requires in 
principle an expression for the r-mode to next to leading 
order in a slow rotation expansion [4]. Following [12, 13] 
we restrict ourselves to the leading order expression for 
the fundamental m = 2 r-mode with density fluctuation 
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where a is a dimensionless amplitude parameter and 
5<I>o denotes the deviation of the gravitational potential 



from the equilibrium. Eq. (4) has a very strong radial de- 
pendence both due to the explicit cubic factor as well due 
to the pronounced r-dependence of the pressure deriva- 
tive of the energy density dp /dp in a neutron star. The 
characteristic time scales of the gravitational instability 
Tq and of the damping by shear and bulk viscosity 
and Tb are defined via their energy dissipation 
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The stability criterium for the r-mode is 1 /xq + 1 /tj + 
1 /%y > where the equality defines the boundary of the 
instability region. At small amplitudes in the subthermal 
regime of the bulk viscosity there is at high frequency an 
instability region where the r-mode is unstable and grows 
exponentially [3, 12, 13]. In the general case the damping 
time due to bulk viscosity reads 
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where / is a constant and the viscosity eq. (2) depends 
on the density fluctuation induced by the r-mode. In the 
integral we neglect the contribution of the crust of the star 
since the bulk viscosity has not been computed there, yet. 

The temperature dependence of the damping times is 
shown in fig. 2 where the bulk viscosity time scale is 
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FIGURE 2. Left panel: Time scales of the fundamental m = 2 r-mode of a 1 .4 Mq neutron star as a function of temperature. The 
dotted curve represents the gravitational time scale, the dashed curve the shear viscosity damping time scale and the solid curves 
the bulk viscosity time scales at various amplitudes. Right panel: Temperature dependence of the "static" saturation amplitude. The 
solid, dashed, and dotted curve are given for the maximum Kepler frequency ~ 6 kHz and lower values. 



given at various amplitudes. As can be seen due to the 
strong increase of the bulk viscosity in the suprather- 
mal regime its damping time decreases strongly with am- 
plitude and eventually undercuts the gravitational time 
scale. Therefore the enhanced damping can effectively 
stop the r-mode growth and saturates the amplitude at fi- 
nite values that are at fixed temperature and frequency 
determined by the above stability criterium. The result 
for this "static" saturation amplitude as a function of tem- 
perature is shown in fig. 2 for different frequencies. It 
features a plateau for temperatures inside the instabil- 
ity region and decreases with frequency. The size of the 
saturation amplitudes should lead to a fast spin down of 
young stars [12]. Actually, during the r-mode growth the 
star cools and starts spinning down so that the above am- 
plitudes do not have to be reached. Due to the continuous 
increase at the right boundary of the instability region 
the r-mode should automatically saturate at an amplitude 
that is sufficient for spinning down the star and that could 
be lower than those of competing mechanisms [8, 9]. 
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